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Управление динамическими системами 
в условиях неопределенности 


Рассматривается решение задач управления, когда качество функционирования динамической систе- 
мы описывается конечномерным вектором показателей качества при информационных предположениях о 
действующих возмущениях на объект и помехах в канале измерения. Предлаг аемый алг оритм управления 
заключается в построении информационного и прогнозирующего множества и нахождении последова- 
тельности управлений из условий совместности системы линейных неравенств. 


Рассматривается решение задач управления, когда качество функционирования 
динамической системы описывается конечномерным вектором показателей качества 
при информационных предположениях о действующих возмущениях на объект и 
помехах в канале измерения типичных для задач гарантированного оценивания и 
управления. Предлагаемый алгоритм управления заключается в построении информа- 
ционного и прогнозирующего множества и нахождении последовательности управлений 
из условий совместности системы линейных неравенств [1-11]. Считается, что в тех 
случаях, когда действия среды перестают быть «враждебными», построена модель 
возмущений [12-19]. В результате уровень неопределенности снижается, что приводит к 
соответствующему изменению множества возмущений. Поэтому в правой части модели 
процессов присутствуют лишь возмущения, известные с точностью до множеств. 
Статья продолжает исследования [1], [8-11], [20]. 


1. Постановка задачи 


Рассматриваются управляемые процессы в традиционной для минимаксного 
[1], [4], [8-10] подхода постановке. Пусть в линейном приближении модель движе- 
ния и измерения задаются системой вида 


ХК! = Ахк + Ви; +Ук, Ук = Чхк1 + Ук, К=0,1,...М- 1, (1) 


где векторы Хк,›Чк› Ук ‚ соответственно, состояния, управления И измерения; векторы 
Хо, Ук, Ук возмущений и помех, известные с точностью до известных множеств 
Хо, М, У ‚, то есть 


ХоЕХу, УкЕУ,, У ЕУ, К=0,1,..., М-1, (2) 


из которых в каждый К-й момент времени возмущения \, и помехи у;.; могут при- 


нимать любое значение. Все матрицы в (1) известны и корректно определены. 
Пусть далее качество функционирования системы (1) оценивается конечномер- 


ным вектором показателей качества ГЕВ“, связанным с вектором состояния ху 
линейным соотношением 


Е = Ехк, (3) 
которое определяет обобщенное множество достижимости (ОМД) [21] 
С: ={ Е =Ехь, хеС, }. (4) 
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Пусть задано множество О; < С; такое, что при ГЕО; функционирование 
системы (1) признается успешным. Поскольку множество достижимости ТФ. полностью 
определяет ОМД С; в К-й момент времени, то далее ограничимся рассмотрением 


построения множеств достижимости. Тогда заданному обобщенному множеству 
достижимости О; будет соответствовать множество достижимости Т.такое, что 


Ух е Г. функционирование системы признается успешным. 
Действительное значение вектора состояния х, является известно лишь в 


виде х; Е Хх, где х,- информационное множество [7-10] 


Хки =Хник ПХ [Ук |, К=0,1..,М-Ь (5) 
Хклик =АХк + Вик + М, Хуа] = [х|бх+у=у,,уеУ}, 


где Хк.лик, Х[Укн | — множество прогнозов и множество совместных с измерением. 
При построении множества прогнозов Х\.!/к выполняется операция суммирова- 

ния множеств АХ; и \, где сумма множеств понимается в смысле Минковского [8] 

(рис. 1), ав (5) производится пересечение множеств Х/к и Х[ Ук. |. Подчеркнем, 


ЧТО ЭТИ действия необходимо выполнять в реальном времени и по существу они 
являются вспомогательными. 


А 


С=А+В 


Рисунок 1 — Сложение множеств в смысле Минковского 


Поскольку значение вектора х, неизвестно, то условие х, ЕТ. заменим на 


условие 
аи (6) 
Множество Г. может также зависеть от К и тогда вместо (6) имеем 
ХС. (к), К = 1,2.....М. (7) 


Условие на множество возможных состояний может быть задано только на 
правом конце 


Хие 74 (М). (8) 
Таким образом, будем выделять три постановки задачи. Необходимо найти после- 
довательность управлений их, К = 1.,... по результатам измерений ук, К= 1.... таких, 


что выполняется одно из включений (6) -— (8). Включение (6) соответствует задаче 
стабилизации, включение (7) — задаче слежения, включение (8) — терминальной задаче. 
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2. Синтез управления 
Рассмотрим решение задачи (7). Обеспечение включения 
Хьник © Фа(К+1) (9) 


будем проводить, представляя множество Та (к+1) многогранником, а множество 


прогнозов Ху.1/к В виде 


0 * 
Хкзик = Хкник + Хкзик + Вх, (10) 
где 
0 * 
Хкик + Хкик =АХк +ГУ. (р 
Здесь р ик — Множество, чебышевский центр которого находится в точке 0, 


* < 
Хк.1/к — Чебышевский центр множества АХ; +ГУ. 


Представим также множество Фа (к + 1) в «центрированном» виде (11) 
(кии = 7. (К+1), (12) 


0 $ * 
где Та — множество, у которого чебышевский центр находится в точке 0, у.| — 
чебышевский центр множества 74 (к + 1) : 


Тогда необходимым условием существования решения задачи (7) является 
выполнение включения 


0 0 
Хьниь +1 © 24 (К+1), (13) 
где Г. — некоторый постоянный вектор. Условию (13) можно придать конструктив- 
ную форму. Пусть у ‚1=0,.. ль, хо, |= 0,...п, — угловые точки множеств г и 


Хки/к Соответственно. Тогда из (13) имеем систему линейных неравенств 


Пт: т: 
о__0 _ ве 
Ув =жм +1; УВЕ =Ь В >0, 1= 1, м, (14) 
т = 
которая должна быть совместной. Пусть условие (13) выполняется для всех К. Рас- 


смотрим случай, когда Г, =0. Тогда из (9), (10) и (12) получим 


* * 
Хкник + Вик = Укнь (15) 
а управление находится соответственно из условия 


* * 
Вик = -Хкнл + кн. (16) 
Пусть теперь Г; #0 и в результате решения последовательности задач (13) 


находится последовательность Г, К = 0,.... Уравнение (15) обеспечивает совпадение 


чебышевских центров множеств Та (к-+1), Хкал. Однако для обеспечения включе- 


ния (9) чебышевский центр множества Х;.1/к необходимо дополнительно сместить 
на Гх. Тогда из (13), (15) получим 


Хкни/к + Вик = Ук — К. (17) 
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Рассмотрим теперь решение задачи (6), которая является частным случаем зада- 
чи (7). Пусть условие (13) выполняется для всех К. Тогда, подставляя в (17) значение 


* < 
чебышевского центра у множества Г ‚ найдем управление в задаче стабилизации (6). 


Для решения терминальной задачи (8) будем строить прогнозирующее мно- 
жество Хи; с 1-го шага на М-Й 


мМ-1 
Хи =Е(М-Ы)Х, + ХЕ(М-Ь}+1(Г\ +В, (18) 
= 
где Е(ь])=А"? 1 <ь Е(Ы) =1 
ИЛИ ВВОДЯ обозначения 
М-1 
Вкал 0; (*) = \Е(М-Ь]+ПВи,, 
= 
М-1 


Гыая \ (+) = > Е(М-Ь]-+П ГУ, 


получим 

Хни=Е(М-М)Х Ги (+) +Вы-ии (*). (19) 

Для центрированных относительно координат множеств х т. . (№) необхо- 
димо проверить следующие условия 

Хи. (М). (20) 

что является необходимым для обеспечения точности терминального критерия (8). 


Если это условие не выполняется, то система оценивания не обеспечивает 
необходимой точности. В этом случае вместо (8) можно использовать критерий 


Хи ПТ, (№) -> тах, (21) 
однако построение управления в этом случае в темпе реального времени не 
представляется возможным, за исключением задач низкой размерности. Компромис- 


* 
сным представляется в этом случае минимизация чебышевского центра Хх 


* 
множества Хх (5) от чебышевского центра уз множества 7. (№) 


хц-7а| > Шщ. (22) 


Пусть условие (20) выполняется и Хх, (№/ 1) +В и: (*) угловые точки мно- 


[к * 


жества Хун (18), тогда из условия совместности систем линейных неравенств 


>. Ву, = х,(№/1)+Вки 9; (+); 
9 (23) 


т: 
УВ, — 1, РВ, >0,]=0,..., т, 
50 


найдем последовательность управлений ц;,....Им—› Доставляющих решение задаче 
Хннс 7. (М). (24) 


Эти управления имеют характер программных управлений. Для управления 
системой используется только первый член и; из найденной последовательности 
управлений, в результате система (1) переходит в состояние хи, которое оцени- 
вается по измерениям у:+1 в виде информационного множества Х;+1. Далее процедура 
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повторяется. Для реализовавшейся позиции хн! строится новое прогнозирующее 
множество Хи! и в результате решения задачи (20) находится последовательность 
программных управлений ц,,..., Иц_| ИТ.Д. 


3. Реализация алгоритма управления 


Предлагаемый алгоритм синтеза терминального управления в условиях неопреде- 
ленности заключается в построении информационного и прогнозирующего множества в 
виде выпуклой оболочки своих вершин и нахождении последовательности управлений из 
условий совместности системы линейных неравенств (20). Причем в ускоренном мас- 
штабе времени необходимо находить вершины множества Хми и решать систему линей- 
ных неравенств (20). Сократить время решения этих задач можно следующим образом. 


Как уже отмечалось, угловые точки центрированного множества ха /; ПОЛНОСТЬЮ 
определяются априорными данными, поэтому еще на стадии проектирования можно 
решить задачу аппроксимации многогранника а д более простым множеством — мно- 


гогранником с меньшим числом вершин, гиперкубом, эллипсоидом. Это позволит 
сократить размерность системы (20). Кроме ТОГО, ИЗ находимой последовательности 
управлений используется всего лишь первый член. Поэтому можно взять неравномерную 


дискретизацию всего временного отрезка Т, = {#1 —&,...„ ща — 2}, АЕЩИ- = 008, 
разбив его на ряд диапазонов Т; = Т. +...+Т; 


2 3 В Каждом из которых временные дискре- 


ТЫ А; постоянны, НО убывают от интервала к интервалу. После построения после- 
довательности программных управлений ц;,....Им На неравномерной сетке на 1-ом 


шаге переходим к расчету на 1+1-м шаге, сдвигая диапазоны «вправо» и отбрасывая 
при этом последний. 

Реализация в реальном времени рассмотренных алгоритмов может предъявлять 
невыполнимые требования к вычислительным ресурсам систем управления 
автономных объектов, в частности ЛА [22]. В этом случае можно воспользоваться 
идеями концепции самоорганизующихся регуляторов [23], [24], что позволяет 
рассматривать модель процесса (1) всего лишь как аппроксимирующую для 
некоторого сравнительно небольшого отрезка времени, что приводит к сущест- 
венному снижению размерности вектора хк и отсюда снижению объема вычислений 
при реализации алгоритмов. Рассмотрим, в связи с изложенным нами, подход к 
операциям над множествами в (5), что приводит к построению иных алгоритмов 
оценивания и управления. 


4. Оценивание и управление и системы линейных 
неравенств 


Пусть по результатам измерений ук решена задача оценивания и построено 


информационное множество, которое может быть задано системой линейных нера- 
венств [8], [9]. Тогда задача оценивания для следующего (К+1)-го момента времени 
может быть сведена к решению следующей системы: 


Оха ТУкы = Ук+в 
ХК -Ахк -Вих -\к = 0; (25) 


Хк ЕХь, Ук Е У, Ук-1 Е У, К = 0, М-1. 
Решая систему (25), находим информационное множество Хки Е ХК. Систему 
(25) можно рассматривать и как неявное задание информационного множества Хки . 
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При развитии эффективных вычислительных методов решения систем нера- 
венств представляется возможным проводить вычисления с накоплением измерений 
и для случаев неизвестных матриц А и В в системе (1) со следующими инфор- 
мационными предположениями 


ХоЕХуо, \кЕ М, У, ЕУ, К=0,М -—1, 


а(,}) [ан (5) ав (63) 55) ®[Ън(Ь), 6 (6) |, 
где Хо, У, У - известные выпуклые компакты, заданные в виде систем линейных 
неравенств; ан (1, ), ав (1), Бн(ь1), Ъв (1,1) — заданные числа. В этом случае из (1), 


(26) 


(25), (26) получим систему нелинейных неравенств вида 


Ок У = У 
Хк+1+1 —Ахки — Вики; — Уи = 0; 


а(5,}) = [ан (Ъ), ав(Ъ) ]; (27) 
(1) ®[Вн(Ь)), в (5) | для всех 1}; 
ХкЕХк, Ук Е У, Ук: Е У, 1=0,..., Г. 


Уравнения фильтра (27) за счет совместной обработки совокупности измерений 
Ук, 1=0,....Г. позволяют не только получить информационное множество Х: для 
момента времени Г, но и решить задачи как сглаживания [8], [9], так и оценивания 
значений реализовавшихся возмущений У; .1.:, Ук, 1=0,...,Г. и параметров а(1.]), 6(.]). 
Если собственные числа матрицы А находятся в круге единичного радиуса, 
тогда начиная с некоторого К > 1+1, 1=0,1,... величинами 
АКх,, АРУ 
можно пренебречь в (4) и аналогично в (6), что приведет к снижению требований к 
вычислительным ресурсам при решении систем неравенств (4) и (6). Это обстоятельство 
имеет принципиальное значение для систем реального времени, поскольку позволяет в 
результате решения системы (6) получить оценки возмущений и использовать их для 
построения адаптивного алгоритма фильтрации. Адаптация позволяет уточнить априорную 
информацию о множествах \\, У возмущений и помех и повысить точность управления. 
Аналогично формулируется задача синтеза позиционного управления [8], [9] 
при наличии фазовых ограничений как задача перевода системы из множества Хкв 
заданное множество Хи в момент времени К + 4, то есть как задача нахождения 
последовательности управлений и. ]=0,..,Ч-—1, являющихся решением системы 


»1=0,1,..,ш-1 при х; ЕХ;, м Е\ 


неравенств вида 


ХК) — АХк+; — Вик; — Мк =0,] = 0,41; (28) 
ХК Е Хк,Хк, Е Хк›Хка Е Хи, Як; Е О, Мк; Е У. 


В (28) в целях упрощения рассмотрен случай, когда матрицы А, В и вектор 
состояния х,,;, ]=0,....4-—1 известны точно. Некоторое снижение требований к вы- 


числительным ресурсам при реализации алгоритмов в реальном времени произойдет, 
если рассматривать ограничения вида 


Укн = Я (%*). Уки = $ (%%), 1=Ьшт, К=0,1,..., (29) 
где ой ЕУ\, У Е\У, К=0,1,... — границы соответствующих множеств \\,У. Известно, 


что любую точку выпуклого множества, не принадлежащую границе, можно представить 
как выпуклую комбинацию граничных точек. Рассматривая выпуклую комбинацию, 
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построенную на нескольких точках, удовлетворяющую условиям (29), можно отказаться 
от выполнения операций над множествами (26). Для этого необходимо построить 
алгоритмы оценивания, идентификации и позиционного управления, работающие только 
с граничными точками вида (29), и за решение задачи оценивания и управления брать 
одну из точек, принадлежащую выпуклой оболочке соответствующего решения. 
Посуществу, это приводит к распараллеливанию алгоритмов, осуществляющих множест- 
венно-множественные отображения на целый ряд существенно более простых алгорит- 
мов оценивания, управления, осуществляющих точечно-точечные отображения. 

Другой подход, позволяющий избежать проведения операций над множествами, 
состоит в том, что возмущения У, Ук,К = 0,1,... в системе (3) можно представить в 


виде выпуклых комбинаций 


Ук = У‘ мк, У а=Ъ а >0, 1=рт; 
11 11 
= У, Вмь УВ =Ь 20, 1=1 м, (30) 
151 11 
точек \', я 1=Ьт, К=0,1,... из (28), где а,б., 1=Ьш в (30) могут рассматри- 


ваться как вероятности, причем \ Е \/, Ук ЕУ, К=0,1,.... Тогда, задавая вероятностный 
механизм генерирования последовательности 

а, =а, (К), В = (К), 1=1, 1, К=0,1,..., (31) 
получим вероятностную аппроксимацию (29) — (31) модели возмущений (26), что 


существенно упрощает реализацию алгоритмов управления динамическими систе- 
мами в условиях неопределенности. 
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В. Ширяев 

Керування динам!чними системами в умовах невизначеностт 

Розглядаеться розв’язання задач керування, коли яксть функщшонування динам!чно! системи 
описуеться ск1иченновимрним вектором показникв якост! при нформашйних передбаченнях щодо 
длючих збурень на об’ект та перешкодах у канал! вимру. Пропонований алгоритм керування 
полягае у побудов! 1нформащйно! та передбачувано1 множини 1 знаходженн! посллдовност! керувань з 
умов сшльно] системи лийних нервностей. 


И’. 1. утауеу 

Оупашес Зу$ет$ Сопы`о! ипдег Соп@ 01$ о? Опсегайиу 

Сопёго]! ргоМет 1$ сопз1Аеге4 г 4упапис зузетз \Возе РЮпсНопте аааШу сап Бе Чезстфе4 Бу а 
Ноице-аптепз1опа| уесюг оЁ диаШу рагатеег$ ипдег поппайопа! аззатрноп$ оп @1атФапсез ап по15е ш 
теазигте зузет. ТВе ргорозе4 сопёго] агогИйт с015155 оЁ зеатсЬ оР пЮптаНопа| ап4 огесазйпе 5е5 
апа Ип4те ой а зеааепсе ог сопёго!$ заб]есеа ю соп4оп$ оЁ сопз1%епсу Юг Фе зу$ет оЁ Ппеаг 
шедиаПиез. 
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